Strong Coupling Effects in Binary Yukawa Systems 
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We analyze the acoustic collective excitations in two- and three-dimensional binary Yukawa sys- 
tems, consisting of two components with different masses. Theoretical analysis reveals a profound 
difference between the weakly and strongly correlated limits: at weak coupling the two components 
interact via the mean field only and the oscillation frequency is governed by the light component. 
In the strongly correlated limit the mode frequency is governed by the combined mass, where the 
heavy component dominates. Computer simulations in the full coupling range extend and confirm 
the theoretical results. 



Recently there has been a great interest in the col- 
lective excitations of Yukawa liquids and solids, cre- 
ated by the emergence of the new field of complex 
(dusty) plasmas |Tj . Complex plasmas consist of highly 
charged mesoscopic grains immersed in the background 
of electrons and ions. It is the presence of the latter 
that, by screening the bare Coulomb interaction between 
the grains, generates an effective interaction that in a 
good approximation can be represented by the Debye- 
Hiickel, or Yukawa potential <f>(r) — Zeexp(—Kr)/r (k 
is the screening parameter). The strength of the cou- 
pling governing the behavior of the systems is conven- 
tionally characterized by the nominal coupling constant 
T = Z 2 e 2 /akT (a is the Wigner-Seitz radius and T is the 
temperature). Due to the screening (k > 0) the effective 
coupling constant T* (defined in [2j, |3j) may be substan- 
tially smaller. The high value of the grain charge (Z>1) 
ensures that the system is in the strong coupling (T* 3> 
1) regime and consequently in the liquid or solid phase. 
Both two-dimensional (2D) and three dimensional (3D) 
Yukawa systems (YS) are of interest, although most of 
the experimental work has focused so far on 2D systems. 

Over the years the collective excitations in YS-s have 
been studied experimentally, by computer simulations 
and theoretically. The analytic treatment of the crys- 
talline solid phase is feasible via the harmonic phonon 
model. More challenging is the appropriate description 
of the liquid phase: here various approximation schemes 
have been attempted [H, HI, out of which the Quasilo- 
calized Charge Approximation (QLCA) Q has emerged 
with considerable success. All these efforts have by now 
congealed in a reasonably complete understanding of the 
collective mode spectra of liquid and solid YS-s, both in 
2D and 3D. In addition, the YS has turned out to be a 
useful paradigm for other strongly coupled many body 
systems 

The restriction in almost all of the foregoing investi- 
gations on collective excitations is that they address one 
component YS-s (OCYS), where all the particles carry 



the same charge number Z and the same mass, m. Only 
a few exceptions are available, most notably the recent 
experimental realization of a two-component (binary) bi- 
layer Q . Thus the collective mode structure of YS-s with 
more than one component, that of the two- component bi- 
nary YS (BYS) in particular, is still an open question. 
It is also by no means a trivial generalization of the sin- 
gle component problem: binary systems are well-known 
to exhibit a wealth of novel physical features: a much 
richer phase diagram, the degree of miscibility of dif- 
ferent phases, new modalities of disorder, the excitation 
of optic modes are among st them (e.g. 0, ED]). The 
connection to problems relating to other liquid or solid 
condensed matter systems are more immediate than in 
the single component case: the listing above provides a 
compass and as systems of expected interest binary ionic 
mixtures, ionic crystals, liquid and solid alloys and semi- 
conductor bilayers come immediately to mind. 

The asymmetry between the two components of a BYS 
is characterized by three parameters: the mass ratio 
mijm\, the charge ratio ZijZ\, and the density ratio 
n-ijnx. In a complex plasma these parameters are not 
independent: most importantly, both the m^/mi and 
the ZijZ\ ratios are determined by the relative grain 
sizes. Theoretically and in simulation models, of course, 
these parameters can be separated; indeed they should 
be so distinguished, in order for one to be able to deter- 
mine the different physical effects brought about by mass, 
charge, etc. asymmetries. We have already shown in 12 1 
that for the purpose of calculating the dispersion rela- 
tion the charge, mass and density ratios can be reduced 
to two asymmetry parameters p 2 = (^2^2) / (Zini) and 
q 2 = (Z 2 m 1 )/(Z 1 m 2 ). 

This Letter addresses the issue of the collective spec- 
trum of a BYS. We study the excitation of the longitudi- 
nal acoustic mode and we investigate how the asymmetry, 
the mass difference in particular, between the two species 
affects the sound speed as the coupling strength T is var- 
ied from the weak coupling T* <C 1 to the strong coupling 
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r* regime. The model is analytically studied both 
for the weakly coupled and strongly coupled liquids and 
for the solid phase: the first is done within the simple 
Random Phase Approximation (RPA) model, while for 
the second we apply the QLCA, whose two-component 
version has already been worked out in earlier works on 
Coulomb systems For the solid phase we assume 

crystalline order and use the harmonic phonon approxi- 
mation. For the liquid phase in the intermediate regime 
between the weak and strong coupling we follow the evo- 
lution of the collective modes as a function of the coupling 
strength through Molecular Dynamics (MD) simulation. 
These approaches yield a consistent picture: in a mixture 
of light and heavy components at weak coupling the two 
components interact through the mean field only and the 
oscillation frequency is governed by the light component. 
With increasing correlations, however, the dominance of 
the light component is diminished: in the strong coupling 
limit the light particles attach themselves to the heavy 
masses and the oscillation frequency is governed by the 
combined mass, where the heavy component dominates. 
As a result, one observes the remarkable and rather dra- 
matic reduction of the sound velocity with increasing cou- 
pling. 

Consider now a BYS, with masses and densities rriA 
and tia (A— 1,2), respectively. Each density can be as- 
sociated with a Wigner-Seitz radius cla and a nominal 
coupling constant Ta = Z\e 2 / ciAkT . The infraction po- 
tentials ar e ip 3 A D (k ) = ^Z 2 A e 2 /(k 2 + k 2 ) and (p 2 A (k) = 
2-nZ\e 2 j\Jk 2 + k 2 . In a OCYS the longitudinal collec- 
tive mode is acoustic at long wavelength (k —> 0), and 
has a coupling dependent sound velocity s 0, [l3| . It is 
the equivalent of this mode in the BYS that we concen- 
trate on. In the calculations we parallel the 2D and 3D 
results. In the weak coupling limit the sound velocity is 
determined by the dispersion relation obtained from the 
RPA dielectric function: 

E/j\0/j \ D/i \ ^j4. 
<PA{k)XA{k,u) ; xa = 2 ■ 
rriA u 

(1) 

For our purpose it is sufficient to use the T=0 "cold fluid" 
result: temperature dependent terms slightly increase the 
sound speed over this value. Eq. (fTJ yields the longitu- 
dinal sound velocity s: 
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uj a B = \J^Z 2 A e 2 nA/m, A , lo 2 a u = ^27rZ A e 2 n A /mAa. 



ujq = y/cjf + w| = yl + P 2 Q 2 is the total plasma fre- 
quency and a—y/aia2. The salient feature of the RPA 
result is that through the total plasma frequency the 
sound speed is governed by the reduced mass. i.e. by 
the light component and, with a sizable difference be- 
tween the light and heavy masses, the presence of the 
heavy component plays a negligible role. 



In the strong coupling limit we calculate the dispersion 
using the QLCA with the input of correlation functions 
obtained from MD simulations. The QLCA is expected 
to provide a reliable description of the strongly coupled 
liquid system, as it is attested by its application both to 
2D and to 3D OCYL-s. (see, e.g. . While we do 

not know how faithfully the QLCA can describe the mode 
structure in binary systems, there is little doubt that it 
is reliable in the long-wavelength [k — > 0) limit. This is 
evidenced by the demonstration [15| that it provides a 
smooth transition to the angle-averaged long-wavelength 
phonon dispersion of a corresponding crystal lattice. 

In the QLCA the central quantity is the longitudinal 
dynamical matrix, C AB (k) dependent on the equilibrium 
pair correlation function Hab{t). We note that since 
masses do not affect equilibrium quantities, with appro- 
priate scaling, all the Hab-s are identical and derivable 
from ho(r), the pair correlation function for the single 
component system. The dispersion relation expressed in 
terms C AB (k) follows from 



\\u 2 5 A B-C L AB {k)\\=Q 



(3) 



We consider first the 3D YBS. In the long-wavelength 
(k — >• 0) limit of interest, the QLCA matrix elements are 



Cfi(k^0) = wf 



(1-^3 + ^ 



Cf 2 (k^0)=u;f 



C&(k-X)) = w? 



(4) 
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Introducing the average charge and mass, 
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Zxm + Z 2 n 2 



mini + m 2 n 2 



Tii +n 2 ni+n 2 
we find that the crucial frequency parameter now is 
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Vp 2 + 1 2 



(5) 



a quantity that has been dubbed in the literature as the 
frequency related to "the pseudo-alloy atom" (FPAA) 
(l6j . or to the average atom in the virtual crystal ap- 
proximation [17 ■HH. From Eqs. (3) and (4) one now 
readily obtains the 3D sound speed as 
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(6) 
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A similar calculation in 2D leads to: 
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e' y h AB (r). (7) 



C/ 3D and £/ 2D in Eqs.(j6]) and (J7J), respectively, are calcu- 
lated from U = [Uu + 2p 2 U 12 + p 4 (7 22 ] /(l + p 2 ) 2 , using 
the corresponding U^g and U^g values. Note that for 
q = 1 we have ui = wn and one recovers a quasi-single 
component behavior [l2j. 

There are two important issues to be noted here. First, 
uj, in contrast to uiq, is governed by the heavier compo- 
nent. Second, the leading terms in the equations for s 3D 
and s 2D are seemingly coupling independent, although 
these equations are valid in the strong coupling limit, 
since the QLCA is a strong coupling approximation. The 
explicit coupling dependence enters through the Uab (r) 
terms only. 

Currently, there is no theoretical understanding of the 
transition bridging the weak, light species dominated, 
and the strong, heavy species dominated, coupling do- 
mains. In order to follow the collective excitations in 
this intermediate coupling region and to verify the pre- 
dictions of the QLCA theory, we generate 3D and 2D 
MD simulation data for the longitudinal current-current 
dynamical structure factors, from which we obtain the 
acoustic speed, over a wide range of T values extend- 
ing well into the crystalline solid region. We analyze the 
liquid phase of several systems, described below, in the 
5 < Ti < T m coupling domain; r m is the coupling value 
where the system solidifies (melts) . The lower limit of T 
is set by the limitations of the simulation technique. 

While the theoretical results are quite general, covering 
the case of a binary Yukawa system for any asymmetry, 
the simulations are confined to systems with mi m 2 
but Z\ = Z%. We study both 2D and 3D systems, both 
with equal densities (n\ = n 2 ) and with unequal densi- 
ties: for the latter we have chosen n 2 = (l/2)ni in 2D 
and n 2 = 3ni in 3D. 

At r > T m , the system is in the solid phase, in 2D 
the underlying lattice structure is hexa gon al; in 3D it is 
bcc or fee, depending on the value of k [20(. With mass- 
asymmetry only, the occupation of the lattice sites by the 
two species would be random, since the potential energy 
is independent of how particles are distributed, and en- 
tropy prefers the disordered distribution. It is only when 
in addition to unequal masses (mi ^ m 2 ), the charges are 
also unequal {Z\ ^ Z 2 ), that a particular crystal struc- 
ture is generated over the lattice in order to minimize 
the energy. Nevertheless, in this work, in anticipation to 
its relevance to realistic systems, we have studied crystal 
structures with maximally symmetric distribution of the 
two components. The corresponding symmetric crystal 
structures are shown in the insets in Fig. 1. (We have also 
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FIG. 1. (color online) Longitudinal sound speed vs. T 1 in 
different systems, obtained from MD simulations. 3D: (a) 
ri2 = 7ii, K = 1, (b) 712/ni = 3, k = 3. 2D: (c) 712 = m, K=l, 
(d) 712/711 = 0.5, k=1. m2/mi=5 for all cases. The dotted 
lines indicate the theoretical low T and high T limits; the con- 
tinuous lines represents the results of the QLCA calculations. 
The panels also illustrate the respective crystal structures in 
the solid phases. 



verified by simulation that the chosen structures would 
be stable configurations for the Yukawa solids, should the 
m 2 /mi ratio be accompanied by Z 2 jZ\ = \/m 2 /mi, as 
it happens in dusty plasmas.) The sound velocities for 
these crystal structures were calculated in the harmonic 
phonon approximation; these were also confirmed by MD 
simulations (of systems with T ~ 10 4 , initiated with lat- 
tice configurations). These values and sound velocities 
obtained from simulations for the disordered phase are 
virtually indistinguishable. 

Our main results for the sound speed are portrayed in 
Fig. 1, where the four panels correspond to the four cases 
listed above. In addition to the results of the MD simu- 
lations, the QLCA prediction for the high V liquid, and 
the results pertaining to the long wavelength phonons 
in the respective crystal structures. For the weak cou- 
pling domain we have indicated a representative value, 
the "cold fluid" sound velocity [Eq.(l)]: inclusion of ther- 
mal effects would only slightly increase this value. Fig. 1 
shows the rather dramatic decrease of the sound velocity 
with increasing T from the weakly coupled to the strongly 
coupled regions. For high T, near r m the agreement of 
the MD results with the QLCA predictions is excellent. 
There is also an almost perfect agreement near T m be- 
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tween these two values pertaining to the liquid and the 
crystal lattice values. The small discrepancy, visible for 
the 3D n 2 — 3nj case can be attributed to the anisotropy 
of the sound velocity in the fee lattice: the liquid results 
correspond to angle-averaged values, while the lattice re- 
sult is given along the chosen {001} direction. This point 
is further elucidated in Fig. 2, where the dependence of 
the sound velocities on the mass ratio is shown for the 
different crystalline solid phases and are compared with 
the QLCA predictions. 

Our calculations show (Fig. 1) that the T-dependence 
of the sound speed is weak and that even at the crystal- 
lization boundary it does not amount to more than 10% 
of the leading term. In view of our earlier statement, it is 
now safe to assume that the validity of the QLCA result 
can be extended to the crystal lattice region. The impor- 
tance of this observation lies in the fact that the notion 
of FPAA has been successfully used for the description 
of phonon dispersion for binary alloys [l6, 17 1 and disor- 
dered systems [lil lijjj. as a heuristic concept. Here we 
have provided an analytic derivation of this behavior. 
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FIG. 2. (color online) Longitudinal sound speed vs. m2/mi, 
obtained from lattice calculations. The labels correspond to 
the cases given in Fig 1. The symbols represent QLCA cal- 
culation results for the strongly coupled liquid phase. 

The physical effect that causes the lowering of the 
sound velocity with increasing correlations may be at- 
tributed to the binding of the lighter particles to the 
heavier ones, resulting in a combined effective mass. This 
latter, defined (e.g. in 2D) by 



mi 



ni +n 2 1 



1 



ni kq (s 2D /ujia) 



(8) 



is displayed in Fig. 3, as a function of the coupling. It 
may be noted that the binding seems to become quite 
pronounced already at the relatively low T=5 value. 

In summary, we have analyzed the behavior of the 
acoustic excitation in a binary Yukawa system consisting 
of two components with different masses, as a function 
of the plasma coupling strength T. Our main focus has 
been to see how correlations affect the way the two masses 
bind into an effective mass. Theoretical analysis at the 
weakly and strongly correlated limits shows that while 
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FIG. 3. (color online) Effective mass vs. Ti, obtained from 
MD simulations, for 2D sytems at different mass ratios and 
densities, at k=1. The dotted lines represent extrapolation 
to low T RPA values. 



in the weakly correlated (r < 1) system the effective 
mass forms as the reduced mass ( "parallel connection" ) , 
in the strongly correlated liquid or solid phase (F > 1) 
the effective mass is simply the weighted average of the 
two masses ("series connection"). As a result, the sound 
speed is substantially diminished in the strong coupling 
domain, as compared to its weak coupling value. Our MD 
simulations of the longitudinal acoustic mode, straddling 
the intermediate coupling range, have confirmed the re- 
sult and have mapped the variation of the sound speed 
over a wide range of T values, both for 3D and 2D sys- 
tems. 

A remarkable feature of the derivation is that the lead- 
ing term in the strong coupling expression is formally 
correlation independent: it is a consequence of the lo- 
calization of the particles, inherent in the model. Thus, 
even though our analysis pertains only to systems inter- 
acting through a Yukawa potential, one may expect that 
it has a more general validity and other binary systems 
with an acoustic or quasi-acoustic type excitation (a two- 
dimensional Coulomb liquid in particular) would follow 
a similar pattern. 

From the experimental point of view, we note that in 
alloys - bearing a great deal of similarity to Yukawa solids 
- the notion of the "pseudo-alloy atom" is heuristically 
well established (see e.g. [l6[ and references therein). As 
to complex plasmas, laboratory experiments on binary 
systems could be more feasible in 3D than in 2D, where 
they should require levitating grains with two different 
Z/m ratios. 
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